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by Kenneth J. Baumelster and Thomas D.Hamill 
Lewis Research Center 

SUMMARY 

The mass evaporation rates and heat transfer coefficients are determined theoreti- 
cally for liquid drops that are supported by their own superheated vapor over a flat hot 
plate. 

The liquid drops are assumed to have a flat disk geometry with a uniform vapor gap 
beneath the drop and a saturated vapor cover. The assumptions are made that the bottom 
of the drop is at the saturation temperature and that evaporation takes place uniformly 
beneath the drop. 

For steady- state laminar incompressible flow, assuming constant properties, the 
Navier-Stokes equations with inertia terms neglected, the continuity equation, and the 
energy equation are solved simultaneously to obtain the heat-transfer coefficient. 

~ * !/4 

\ = 0. 68(k J A*gp £ p v /ATfxL e ) 


For heat transfer to the drop 


A* = A[l + (7C p AT/20A)] ” 3 


For heat transfer from the plate 

A* = A [l + (7C p AT/20A)] 5 J [l + (l/5)(C p AT/A)] 4 

where h is the heat-transfer coefficient, k is the thermal conductivity of the vapor, 
p is the density, {i is the absolute viscosity, L is the equivalent geometry factor, 
and X is the heat of vaporization. 

These closed form results agree with the value of the heat-transfer coefficient pre- 
viously obtained by a computer solution of the exact Navier-Stokes equations. The earlier 
work was in excellent agreement with experiments. The conclusion is drawn that inertia 
terms in the Navier-Stokes equation are negligible relative to the viscous terms. 



INTRODUCTION 


Investigations of heat transfer to liquid drops that are supported by their own super- 
heated vapor were begun as far back as 1756 when Leidenfrost (ref. 1) described the 
phenomenon of film boiling. At present, references 2 and 3 indicate considerable interest 
in the subject of drop vaporization in the fields of internal combustion engines and metal- 
lurgy. 

Dimensional and semiempirical correlations for the evaporation of liquid drops in 
film boiling on a flat plate are presented in references 4 and 5. Reference 6 presents a 
solution of the exact Navier-Stokes equation and the energy equation; however, a grapi- 
cal approach was used to combine the energy and momentum equations. Although this 
approach gives a physical insight into the coupling between the momentum and energy 
equations, the relation of the heat-transfer coefficient to the basic fluid properties is 
obscured. 

Reference 7 shows how the computer solution to the exact Navier-Stokes equation 
and energy equation can be combined to yield a closed form solution for the heat- 
transfer coefficient h. The resulting form of the heat-transfer coefficient is similar in 
nature to that derived by Bromley (ref. 8) for heat transfer in stable film boiling. 
Bromley’s analysis is based on the solution of the momentum equations in which the 
inertia terms have been assumed negligible. 

The similarity between the results of references 7 and 8 led to the observation that 
a closed-form solution of the governing equations without the use of computer solutions 
could be obtained by neglecting inertia effects (creeping flow). 

Solution of the Navier-Stokes equations neglecting inertia terms is obtained, and the 
results are compared to the computer solutions, which had considered the inertia effects. 


SYMBOLS 

A area, ft 

a constant of proportionality, sec”^ 

C constant 

0^ specific heat of vapor at constant pressure, Btu/(lb mass)(°F) 

e z unit vector in upward z-direction 

F force, lb force 

f transformation variable, ft/sec 
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axial pressure variable, sq ft 

o 

acceleration of gravity, ft/sec 

dimensional conversion factor, 32.1739 (^)( lt5 mass )/(l b f orce )( s ec^) 

heat-transfer coefficient, Btu/(hr)(ft^)(°F) 

thermal conductivity of vapor, Btu/(hr)(ft)(°F) 

equivalent geometry factor, ft (see eq. (31)) 

average calculated drop thickness, ft 

mass of drop, lb mass 

2 

static pressure, lb force/ft 

o 

environmental pressure, lb force/ft 
Prandtl number 
rate of heat flow, Btu/ (hr)(ft^) 
radius, ft 

maximum radius of drop, ft 
temperature, °F 



time, sec 

radial velocity, ft/sec 
drop volume, ft 
velocity, ft/sec 
axial velocity, ft/sec 
axial coordinate, ft 

2 

thermal diffusivity of vapor, a = k/p C , ft /sec 

v p 

gap thickness, ft 
z/ 6 

heat of vaporization, Btu/lb mass 
modified heat of vaporization, Btu/lb mass 
absolute viscosity of vapor, lb mass/(ft)(sec) 
kinematic viscosity of vapor, ft /sec 
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£ dummy variable, ft 

Q 

p density, lb mass/ft 

r dummy variable, ft 

Subscripts: 
f film 

i index 

£ liquid 

p plate 

sat evaluated at saturation condition 
v vapor 

z axial coordinate 

6 evaluated at lower surface of drop 

Superscripts: 

* derivative with respect to independent variable 
— vector 


BASIC MODEL AND EQUATIONS 

The model used in this study applies to large drops for which a flat disk model 
(fig. 1), can be used. For example, references 5 to 7 indicate that for water drops in 
the volume range 0. 05 to 1 cubic centimeter, the analytical model based on a flat disk 
geometry reasonably satisfies the physical situation. 

The following assumptions are made in developing the analytical model: 

(1) Heat transfer and evaporation from the upper surface are considered negligible 
in comparison to rates of transport beneath the drop. 

(2) Radiation is neglected. 

(3) Physical observation indicates that the gap thickness is a function of radial posi- 
tion; however, to make the problem mathematically tractable, a uniform gap thickness 

6 is assumed. 

(4) In a similar manner, the geometric shape of the drop is approximated by a flat 
disk with constant thickness l, where It is defined by the equation 
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Figure 1. - Schematic model of evaporation of flat disk. 


The relation between V and r Q 


v = V(r Q ) 


( 2 ) 


was found numerically by balancing the gravitational force against the surface tension 
force (ref. 6). 

(5) Because of a low calculated Reynolds number (ref. 6), the flow is assumed lami- 
ner and incompressible with negligible energy dissipation. In addition, the body force of 
the vapor in the momentum equation has been neglected. 

(6) The inertia terms in the Navier-Stokes equation are neglected. Justification 
for this assumption is given in the main text of the report, as well as in appendix A. 

(7) In the energy equation, it is assumed that 


u^i« w 
3r 


3T. 

3z ’ 


!_i(r 9 l)« 9 2 T 

r 0r \ 3r/ -2 

' ' az 


( 3 ) 


(8) The velocity and temperature profiles at any instant are assumed to be in steady 
state. 

(9) The bottom of the drop is at the saturation temperature. 

(10) Evaporation takes place uniformly beneath the drop. 

(11) The properties of the flow field are evaluated at the film temperature, defined 


I 


5 


(4) 



+ T 


sat 


2 


Consequently, for this case of axisymmetric flow, the governing equations are as 
follows: 


(1) Momentum, 


„ g c ap 
0 = + v 

P v 9r 


<d\ i3u u_ afu 

. 2 + r ar 2 .2, 

,3r r 3z 


n _ g c ap ^ / a 2 w x law, a 2 w\ 

U — — + V I + *r 


p 3z 
H v 


^3r 


r 3r 


3z 


(2) Continuity, 


(3) Energy, 


*i + H + d JL = 0 

3r r 3z 


for the boundary conditions 


w 3 -I=a^l 


3z 


3z 


z=0 u=0 w=0 T=T 


P 


z = 6 


u = 0 


w = w(6) 


T = T 


sat 


r = 0 u = 0 


r=r Q z - 0 P = P Q 

(4) Static-force balance (neglecting reactive force, see appendix B), 


(5) 

( 6 ) 


(7) 


( 8 ) 


(9) 

( 10 ) 

( 11 ) 

( 12 ) 
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( 13 ) 



P(r, 6)2nr dr = Vp £ 3 - + 7rr^P 

O n 


(5) Interface energy balance, 


- P y Xw( 6) = 



(14) 


The static-force balance has a physical interpretation; in that the weight of the drop 
is assumed equal to the force resulting from the pressure difference between the upper 
and lower surface of the drop. Mathematically, equations (13) and (14) are necessary to 
make the problem determinate, since the values of 6 and w(6) in boundary condition 
(eq. (10)) are unknown. 

The solution involving the above set of differential equations and boundary conditions 
will now be performed. Figure 2 presents a solution flow chart of the problem. This 
chart may be used as an aid in understanding the analysis now presented. 

The differential equations (5) to (8) contain two independent and four dependent vari- 
ables. The boundary conditions specified by equations (9) to (14) do not make the problem 
determinate at this time. For example, additional boundary conditions on w are re- 
quired at the radial boundaries; however, as seen in figure 2, a similarity analysis is 
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Figure 2. - Solution flow chart 
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used in the solution. The similarity transforms will reduce the order of the equation and 
will satisfy the boundary conditions so that the problem becomes determinant. 


METHOD OF SOLUTION 
Momentum Equations 

By assuming constant fluid properties, the interaction of the equations of motion 
with the energy equation ceases, and the velocity field no longer depends on temperature. 
A similarity transform is now used that reduces the partial differential equations 
(5) to (7) into a set of ordinary differential equations. Since the governing equations and 
boundary conditions between the three-dimensional stagnation flow problem and this 
problem are similar, the transforms for the three-dimensional stagnation flow problem 
are used (ref. 9): 


w = -2f(z) 

(15) 

u - rf’(z) 

(16) 

2 

P - 1 Pv Ir 2 - r 2 + G(z) 

2 g r L° J 

(17) 


These transformations satisfy the requirements of reducing the partial differential 
equations in two variables into an ordinary equation in one variable and of satisfying the 
radial boundary condition of equation (11). In addition, these transformations reduce the 
functional dependency of the dependent variables on the independent variables in such a 
manner that the problem is now determinant. For example, w is no longer a function of 
r and thereby requires no specification at the radial boundaries. 

Substituting equations (15) and (16) into equations (5) to (7) results in the ordinary 
differential equations 


f»" = _fL_ (18) 

v 


and 

G* = — f<* (19) 

2 

a 
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Solving equation (18) and applying boundary conditions (9) and (10) result in 


f = 


6 v \2 



( 20 ) 


But expressing f in terms of w by equation (15) yields 


w = -! — (— 6z 2 - z 2 ) 
!i\l / 


Therefore, at the bottom surface of the drop 


w(S) = - liL e 3 
6 v 


( 21 ) 


( 22 ) 


Thus, the axial velocity at the surface of the drop is directly proportional to the cube 
of the gap thickness. Physically, for a given pressure distribution under the drop less 
flow is required to keep the same pressure level for smaller values of 6, as prescribed 
by equation (15). The constant a is now determined from conditions of static equilib- 
rium. 

Integrating equation (19) and substituting in the first derivative from equation (20) 
result in 


G = -2(6z - z 2 ) + C 

Substituting into equation (17) gives the pressure distribution as 


(23) 


p 4 a2 r : [ r “' r2 ' 2 ( 6z '^ + c ] 


(24) 


Solving for C using boundary condition (12) and substituting the value of C into 
equation (24) result in 


1 ? P 
P = - a 2 


vf 

r o 



(25) 
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At the bottom surface of the drop 


P(r,6) = |a 2 f!:(r2.r 2 ) + P 0 

Op 


(26) 


Substituting the above relation into equation (13) and solving for a z result in 


a 2 = i v/— 


77 Wr 4 


(27) 


Therefore, the final expression for pressure is 


2 Vp j l 


P(r, z) = P +- 

77 r 4 ^0 


A(r 2 -r 2 )-i^JLfe Z -z 2 ) 
g.A ' n _4 sr ' ' 


r o ° c 


(28) 


When equation (27) is substituted into equation (22), the axial velocity at the surface 
becomes 


w(6) = -7IV 


2 § 3 


Pv V 


(29) 


Note that for the flat disk geometry used in this model 


V = TiryL 


(30) 


Defining 


V A 2 


yields 


(31) 


w(6) = -AZJL _g_ 6 3 


3 ffp y vL, q 


(32) 
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Thus, for the drop to remain in equilibrium, the magnitude of the axial velocity 
leaving the surface of the drop is proportional to the cube of gap thickness. 


Energy Equation 

Since w is a known function of z, the energy equation can now be solved directly. 
Equation (8) can be rewritten in the form 


A In® = 3 
dz dz a 


(33) 


Integrating the above equation and applying the boundary condition at z = 0 (eq. (9)) 
esult in 


Z d£ exp /^^Ildr 

Jo a 

At z = 6, 

d| exp dr 

Jo 01 

Therefore, dividing equation (34) by equation (35) yields 



Substituting the value of w from equation (21) into equation (36) gives 
integration 




(34) 


(35) 


(36) 


for the first 


exp 



^HdT 


= exp - 



Ova 




(37) 
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Since £ < 6 


/ W < T > dT 

< 

2*4 
a 6 

/ 01 


12^a 

«4) 




( 38 ) 


As a general numerical example, for a 1- cubic centimeter water drop on a flat plate 
at 600° F (appendix A, and noting that a = v/Vr) 


aV 

12va 


= 0. 056 


(39) 


Therefore, the exponential may be expanded and the second and higher order terms 
neglected. Thus, 


exp 



^lidT 2 1 - 

a 



(40) 


Substituting equation (40) into equation (36) and continuing the integration result in 


T - 
T sat 





av 40 / 


(41) 


The temperature gradient of any point z is found by a simple differentiation of 
equation (41) resulting in 


dT 

dz 


^sat 



1 



iz 4 - 

\60 


J_a^ 
40 av 



(42) 


At z = 6 


12 



(43) 


or at z = 0 


dT 

dz 


. < T sat - V 


z=5 


a 2 4 

1 - J* 6 * 

12 av 

1 _ JL 6 4 

40 av 


dT 

dz 


T sat ~ T p 


z=0 


i-i. Sle 4 


40 av 


(44) 


The above expression for the temperature gradient will be used in the determination 
of the heat-transfer coefficients in the next section. 


Heat-Transfer Coefficient To The Drop 


An energy balance on the liquid drop results in 


-p v w(6)X = 



z=6 


(45) 


which is the interface energy balance of equation (14). 

Substituting the value of w from equation (22) and the temperature gradient from 
equation (43) results in 



(46) 


When it is noted that 
theorem gives 


2 fi 4 
a 5 

12 ^ 0 ! 


is small compared with 1, expansion by the binominal 
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( 47 ) 


2 /i 

1 -^— 5 4 
12av 

40 ctv 


1 -JL£ 8 4 

120 ay 


as an approximation when second and higher order terms are neglected. 
Solving equation (46) for 6 results in 




< T p - T sat> k 


Substituting the value of a z from equation (27) yields 


r 


i 

p< 

.H ! 

a 2 X 

rt 

r 

1 + - 

a 

^ v 6 

2 


~\ 


sat' 


1/4 


=r > 


6 = J 


r 


3 ’ k < T p - T S at>" r o 

2p v p £ VgX 

V 

r, . 7 C p (T P ' T sat> 

20 X 


1/4 


The heat-transfer coefficient to the drop is defined as 

q ( - -P/w<a) = h t (T p - T sat ) = h f AT 

Therefore, 

-p v Xw(6) 


h. = 


AT 


(48) 


(49) 


(50) 


(51) 


Substituting the value of w from equation (22), the value of 6 from equation (49) 
and rearranging terms result in a heat-transfer coefficient to the drop of the form 

*3 ,>* XV4 

k P v p f gX £ 


h £ = 0. 68 | 


fi ATL 


(52) 
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where the modified latent heat of vaporization X^ is defined as 


>* - 


( 53 ) 


7 C n AT 

|1 + _L 

20 X 


Heat-Transfer Coefficient From Plate 


The heat-transfer coefficient from the plate is defined by the following equation 


h AT = -k — 
p 9z 


(54) 


z=0 


Substituting in the value of the temperature gradient given by equation (44), using equa- 
tion (49) in the first term of equation (44), and rearranging terms give 


h = 
P 


k 3 p v a 2 X (1 + 


7 C P AT 
20 X j 


6 ATy 


1/4 


' 9 ' 

1 - — — 5 4 


(55) 


40 av 


Substituting the values of a (eq. (27)) and L (eq. (31)) into the above equation yields 


/ k 3 p p £ gX 

h = o.68 — 
p \ p ATL 


1/4 


P 


(56) 


where 


A 1 + 


* ' 
A* = - 


7 C p AT 
20 X / 


l-i^6 4 f 

. 40 av / 


(57) 
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But substituting into equation (57) the value of 6 from equation (49), and the value of a 
from equation (27), and clearing fractions yield 


X 


* 

P 



+ 


_7_ 

20 



CpATj 4 


(58) 


or 


X 


* 

P 



7 C p AT \ 5 
20 X j 


1 

+ — 
5 




4 


(59) 


The difference between and results from sensible heating effects in the 
superheated vapor. 

Equations (52) and (56) agree with the functional form of Bromley T s equation (ref. 8) 
for film boiling heat transfer from a horizontal tube. The only significant difference is 
that the diameter of the tube is replaced by the newly defined geometry factor L 0 . 


COMPARISON TO COMPUTER SOLUTION OF THE EXACT NAVIER-STOKES EQUATION 

References 6 and 7 present solutions to this problem by solving the exact Navier- 
Stokes equation. It is enlightening at this time to compare some of the results of this 
earlier work, which considered inertia effects, to the results of the present paper, which 
has neglected inertia effects. 

It was observed in reference 6 that the radial velocity profile was for all practical 
purposes parabolic in shape. This agrees with the result of the present paper, since 
combining equation (16) and the first derivative of equation (20) yields 

u = - - — - (z 2 - 6 z) (60) 

2 v 
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or 


u = 


ra 2 5 2 / 2z 
8v \6 



(61) 


which for a fixed r, indicates that the velocity distribution is parabolic in shape. 

The axial velocity at the surface of the drop, as given in reference 8, is 

a 2 8 

w = -0. 172 — — 5 3 (62) 

v 


Likewise, this compares very closely with the results shown in equation (22). In fact, 
the difference in these two coefficients could be due mostly to computer inaccuracies 
rather than to the magnitude of the inertia terms: except for corrections on the modified 
latent heat of vaporization in equation (53) and (59), the expressions for the heat-transfer 
coefficients are equal to the results in reference 7 where inertia terms were considered. 

Finally, a numerical calculation is presented in appendix A which shows that the 
radial and axial components of the inertia terms are negligibly small at the surface of 
the drop and at the plate; however, the inertia term for the axial component is significant 
in the center of the vapor gap because the viscous term goes to zero. Nevertheless, 
conditions at the center of the gap do not control the heat-transfer process going on at 
the wall and at the drop surface. Consequently, the inertia effects can be dropped in 
this analysis. 


CONCLUDING REMARKS 

In the evaporation process to a liquid drop, the heat-transfer coefficient is shown to 
be equal to 


h. = 0. 68 


3 * \ 1/4 

fc 3 A*gp £ p v ) 

AT pL 


For heat-transfer to the drop 


* * / 7 C p AT 

1 \ 20 A , 


17 



For heat-transfer from the plate 





~ 


5 




4 


It is shown by a comparison with a numerical solution of the exact Navier-Stokes 
equations, that the inertia effects in the flow field are negligible. The analytic expres- 
sions for h^ are identical in form to the Bromley equation for film boiling from a 
vertical plate. Theoretical and experimental values of the heat-transfer rates are in 
good agreement as shown in reference 7. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, August 10, 1965 
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APPENDIX A 


COMPARISON OF INERTIA TERMS WITH PRESSURE AND 
VISCOUS TERMS IN NAVIER-STOKES EQUATION 

The validity of the creeping motion approximation can be checked by evaluating the 
acceleration terms in the momentum equations using the velocity functions (eqs. (15) 
and 16)). 

If the acceleration terms are small compared with viscous or pressure terms, at 
the very least, the solutions are consistant with the assumptions under which they were 
derived. 

Inertia effects in the radial direction are first considered. The ratio of viscous to 
inertia acceleration in the radial direction is given by 


-2 
3 u 


3z 


.2 r 


u?H + w^ 
3r 3z 


2*4 
a o 


1 , 2\^" l/i o \/3 2 3 

- (V ~ V ) - - (1 - 277) -?7 - t] 

4 6 \2 


(Al) 


As a typical case, consider a 1- cubic- centimeter water drop on a plate at 600° F and 
atmospheric pressure. For this case (ref. 6) 


a 2 = 4. 5xl0 6 


sec 


-2 


v = 4. 0x10" 


ft 2 /sec 


6-4. 0x10' 


ft 


By using the preceding equations, the ratio expressed in equation (Al) has a minimum 
value at 77 = 1 of about 18, and goes to infinity at the wall. Thus, the viscous terms 
always dominate the inertia terms for the radial flow. A plot of the ratio (eq. (Al)) is 
shown in figure 3. 

Finally, consider the axial component of momentum with inertia terms retained 


w dw = _ g cdV + ^ 

dz p ?z dz 2 


(A2) 
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Figure 3. - Comparison of viscous with inertia terms in 
radial direction. 


where w = w(z) only. 

It must be shown that 




dz 


dz 


(A3) 


Using the expression for w(z) 

(eq. (21)) and defining rj = z/5 yield 

dw a 4 6 5 /3 2 3V 2\ 

w — = -77 - rj Jl 77 - 77 J (A4) 

dz 3^2 V 2 A ) 

and 



a 2 6(l - 2rj) 


(A5) 


For the 1- cubic centimeter drop 
previously considered the preconstants 
of the 77 functions become 


a 4 5 5 = 430 ft/sec 2 


a 2 6 = 1800 ft/sec 2 


(A6) 


(A7) 


The inertial acceleration in the z-direction w dw/dz goes through a maximum at 
about 77 = 0. 686. The value of the inertial acceleration is computed to be 35.3 feet per 
second squared. At the same distance from the wall the viscous term is about 20 times 
greater than the inertia term; however, the results of this calculation can be misleading 
because there is a region near the center of the gap where the computed inertia terms are 
larger than the viscous terms. For example, at z = 6/2 the viscous acceleration term 
is identically zero; whereas, the inertia term w dw/dz is equal to 27 feet per second 
squared. A complete picture of the region where the inertia terms become appreciable 
is given in figure 4, where the ratio of viscous to inertia terms is plotted as a function 
of the distance from the plate. In the range 0.44 <z /6 < 0. 60 the ratio of viscous to 
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Fractional distance from plate, z/6 


Figure 4. - Comparison of viscous with inertia terms in 
axial direction. 


inertia terms becomes less than 10 in 
absolute magnitude. Outside this region 
the inertia forces are small compared 
with viscous forces. 

The region of inconsistency encoun- 
tered in comparing viscous and accelera- 
tion terms in the axial direction might 
cast some doubt on the analytical solu- 
tion, Ixit it is felt that the only effect of 
the inertia terms is to distort the velocity 
profile somewhat in the center of the 
vapor film. 

The inertia effects on the velocity, 
pressure, and temperature profiles in 
the vicinity of the wall or the vapor- 
liquid interface are small. The heat- 
transfer coefficient depends mostly on the 
thickness of the gap, which is determined 
by the pressure distribution at the liquid- 
vapor interface. Since the inertia effects 
produce only minor perturbations of the 
pressure field in the vicinity of the inter- 
face, the computed gap thickness and 


heat-transfer coefficient should not be affected by neglecting the inertia terms. In the 


problem of film boiling in a forced- convection boundary layer, a similar conclusion was 


reached in reference 10: inertia forces are not important in the vapor film for the range 


of parameters usually encountered in practice. 
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APPENDIX B 


REACTIVE FORCE 

In reference 11, the momentum theorem for a control volume V is given by 

Z ? = * ///? dV + ft % m 

All the elements of the drop were assumed to move at the same velocity 

J JJ g c g c 

For this problem and W y are parallel to the z-axis and may be conveniently 
expressed as the product of a scalar and a unit vector e z in the axial direction. Thus, 


w v = w v (6)e z 

(B3) 


(B4) 

dA = -dAe z 

(B5) 


Therefore, the force balance in the axial direction yields 


F 


z 


3M . VT £v w 2 (5)dA 

Sc 51 S c 8t W Sc 


(B6) 


If the gap thickness is assumed to remain approximately constant over a short period 
of time, and W y can be considered constants. Thus, 



g c 3t 


- ^ W y (5)A 
g c 


(B7) 
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(B8) 


dM 

dt 


= P V AW V (6) 


Therefore, 


£ F ,=^K- W v (6 >] (B9) 

If it is assumed that the drop moves uniformly downward to keep the gap thickness 
constant, 


dM 

W„ = (BIO) 

* p 4 a 

Therefore, substituting equations (B8) and (BIO) into equation (B9) results in 

For regions sufficiently far from the critical point, the density ratio is negligible 
compared with 1; thus, 



PyAwfo) 


(B12) 


The forces exerted by the surroundings are a result of pressure and gravity. Thus 


-VP, f + 


r 


P(r, 5)27rr dr 


" r o p o 


Pv*w|(6) 


(B13) 


The term on the right side of equation (B13) is defined as the reactive force assoc- 
iated with the ejection of mass from the drop surface. 



For a plate temperature of 600° F and a 0. 5- cubic centimeter water drop, the re- 
* 

active force is computed (ref. 6) to be 4. 6X10 pound force; compared with the weight 
of the drop, 1060x10” pound force, the reactive force is not important. Thus. 



P(r, 6)277r dr = Vp £ 


t * ’ r ° p ° 

®c 


(B14) 


Equation (B14) is referred to as the static-force balance in the text (eq. (13)). 
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